The effects of spin-spin interactions on magnetoresistance in disordered organic semiconductors 
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A recent theory of magnetoresistance in positionally disordered organic semiconductors is extended to in- 
clude exchange and dipolar couplings between polarons. Analytic results are discovered when the hyperfine, 
exchange, and dipolar interactions have little time to operate between hopping events. We find an angle-of-field 
dependence of the magnetoresistance that agrees with previous experiments and numerical simulations. In ad- 
dition we report new magnetoresistive behavior that critically depends upon the amount of anisotropy in the 
dipolar interaction. 



I. INTRODUCTION 

Some spintronic applications are hindered in many inor- 
ganic semiconductors due to the heavy constituent atoms 
present and the resulting strong spin-orbit coupling that 
severely limits spin lifetimes. The less massive constituent 
atoms of some organic systems offer the advantage of weak 
spin-orbit coupling as well as other benefits such as cheap 
processing and chemical tuning. ' ^ In recent years, spurred by 
magnetoresistance measurements of 10-20% in small fields at 
room temperature, the effects of magnetic fields on currents 
through organic layers have been studied by many researchers 
with hopes of adding magnetic functionality to existing and 
new organic devices.^ '^ This large magnetoresistive phe- 
nomenon, known as organic magnetoresistance (OMAR), has 
typically been explained in terms of spin correlations between 
charge caiiiers or polaron pairs. There is strong evidence^^ 
that OMAR occurs whether the spin correlations are between 
like charge carriers"^ '^ '^ or oppositely charged carriers.^ ^ 
By taking note of the intimate relat ionsh ip between spin and 
charge transport,'^ two recent work a'^ l '^ ^have proposed a the- 
oretical description of OMAR based on percolation theory in 
positionally disordered unipolar organic systems. The goal of 
the current article is to apply this theory to include exchange 
and dipolar interactions between polaron pairs and thereby es- 
tablish the wider applicability of the theory. 

Wagemans et al. studied the problem of incorporating spin- 
spin interactions experimentally, by measuring the field an- 
gle dependence of OMAR, and theoretically, by using the 
stochastic Liouville equation approach.'** Prior to the mea- 
surements of Wagemans et al. no field angle dependence in 
OMAR had been reported, though in 1970's magnetic field 
angular dependence had been observed in the photolumines- 
cence of an organic crystal. '^ Wagemans et al. also posit that 
anisotropic hyperfine interactions can lead to similar angular 
dependence in OMAR behavior. Anisotropic hyperfine inter- 
actions are not considered in this article. 

We apply spin-dependent percolation theory to organic 
magnetoresistance in the presence of exchange and dipolar in- 
teractions and (1) obtain analytic expressions for the angular 
dependence of the OMAR, including the number of extrema 
as a function of angle, in the regime of fast hopping, (2) pro- 
vide a heuristic picture of the effect of exchange on OMAR, 
by inferring that it behaves as an effective hopping for spin 
pairs, and finding an analytic correspondence between hop- 



ping and exchange valid for all regimes of hopping, and (3) 
predict large negative MR for special angles where the spin 
splitting from exchange and dipolar interactions vanishes in 
the high-field regime. In the regime where analytic expres- 
sions can be obtained, we find in the limit of hopping times 
going to zero that the angular dependence vanishes (in agree- 
ment with the numerical results of Ref. [TSl l. This theory also 
explains, via analytic expressions, how two unique angular 
dependences of the saturated magnetoresistance develop for 
certain ratios of the exchange and dipolar constants, and de- 
scribes the interplay between them when both are present; in 
all three scenaiios magnetoresistive behavior distinct from the 
non-interacting case is seen. Notably, spin-spin interactions 
allow for the presence of an ultra-small-field-effect (USFE), 
where the magnetoresistance changes sign in small fields. We 
further find that the amount of anisotropy in the dipolar dis- 
placement vector is crucial in determining the angle depen- 
dence of the magnetoresistance. 

Section |ll] describes the semiclassical treatment of the hy- 
perfine fields along with the exchange and dipolar spin-spin 
interactions, and for completeness of this article Section III 
summarizes the percolation-based model of OMAR from 
Refs. [T5]and[T6l This model is then applied in Section [rV| to 
the fast hopping regime of caiTier transport in organics in the 
presence of exchange and dipolar spin-spin interactions, and 
analytic results are derived. Due to the presence of spin-spin 
interactions, the theory outlined in Ref. 16 for OMAR in a 
regime of arbitrary hopping times must be extended; sectionM 
introduces a different formaUsm that can include spin-spin in- 
teractions and arbitrary hopping rates, and then describes the 
full magnetoresistance problem with arbitrary hopping rates 
and inclusion of exchange and dipolar spin-spin interactions. 



II. THE HYPERFINE, EXCHANGE, AND DIPOLAR 
INTERACTIONS 



Consider a polaron pair (PP); the physical picture is two 
spin- 5 polarons located at two sites. The spins evolve co- 
herently under the influence of identical applied fields, differ- 
ent classical nuclear or hyperfine fields, and the exchange and 
dipolar fields. The respective PP Hamiltonian is 






(1) 



whose terms are expressed later in this section. All Hamil- 
tonians and spin operators are written in units of h through- 
out. Since nuclear dynamics are slow, the total hyperfine field 
at a site is stationary and composed of many nuclei (Figure 
lljl, which justifies the classical approximation.^*' The hyper- 
fine fields are random site-to-site. When a polaron hops to 
an unoccupied site (disassociation), its coherent spin evolu- 
tion ceases as the polaron now feels a new local magnetic 
field. When the polarons hop to combine as a bipolaron, they 
are necessarily in the singlet state and the large exchange in- 
teraction prevents further spin evolution. In the PP, the two 
polarons are further apart, exchange is much smaller which 
allows spin transitions. If the hopping is fast, the PP spin 
evolves very little due to the motional narrowing effect. When 
the spins reside on their sites for a long period of time, they 
precess around their hyperfine fields many times and are "fully 
mixed". 
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FIG. 1. Semiclassical description of total field. Each site is com- 
posed of some number of spin-moment-carrying nuclei which com- 
bine with the external field (Oq to create ©ai. Each site has a different 
total field picked from the distribution of Eq. ([3}. 



A. The hyperfine interaction 

The Larmor frequency due to the conglomerate nuclear spin 
is the constant classical vector 
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In Section III we discuss quantities that require hyperfine 
averaging. Explicit averaging by integrating over the six hy- 
perfine variables is a time-consuming numerical task. We opt 
to average by considering many different hyperfine configura- 
tions taken from the Gaussian distribution. Typically we use 
10^ — 10^ configurations. 



where Oj is the hyperfine coupling constant between the elec- 
tron and the y-th nucleus. I^., the precession rate for carrier 
at site ;, is made up of many different nuclei each with vec- 
tor length a j ^yJj{Ij~+Y) pointing in a random direction. /, 
is a spin quantum number (1/2 for a proton). The probabil- 
ity distribution for finding an / 'molecule' among an ensem- 
ble of such molecules with its total end-to-end vector between 
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is the effective hyperfine coupling width due to all the nuclei 
at a site. We consider all sites to be of the same type. The 
effective ^eW is Be/f = Oeff/je with je = 0.176 ns^' mT^^. 
Beff is on the order of 1 mT. 

A polaron spin precesses at a frequency 



WtV,- =/iV,+(Oo, 



(5) 



where ©o = JeBz is the applied field. When neglecting spin- 
spin interactions, the non-interacting Hamiltonian of interest 

is 

^Q = ^4/ + ^^z = {In, + ©o) • 5i + {In2 + ©o) • ^2 (6) 
whose terms can be identified in Figure [T] 



B. The exchange and dipolar interactions 

Due to the Coulomb interaction and the quantum mechan- 
ical requirement of an antisymmetric total wavefunction, two 
polarons experience a mutual exchange interaction: 



Jifixc^J{r)d+2Si-S2), 



(7) 



where 2J{R) is the exchange coupling constant which defines 
the T-S energy splitting (in frequency units). J{R) decreases 
exponentially with distance between the spin pair. Since PPs 
have a large assortment of separations (R), the exchange inter- 
action is expected to vary widely throughout a array of PPs. 

Spin carrying polarons in the vicinity of one another expe- 
rience each other's spin magnetic moment. For simplicity we 
assume point dipoles. The energy of the interaction depends 
on the orientation of the spins. This magnetic interaction is 
known as the dipolar interaction with the Hamiltonian 



Jfdip^D{Si-S2-3iSi-R)iS2-R)) 



(8) 



where D = iAQg^ii\/A%hP? . R is the vector displacement 
between the two polaron spins. If R was uniformly dis- 
tributed between PPs, one does not expect any OMAR an- 
gular dependence;^^ however there is reason to believe that 
there is an anisotropy in the formation of PPs such that R is 
restricted to a narrower range of directions; specifically R may 
align with the sample normal as shown in Figure Etc). ^^ This 
is a sensible expectation since an electric field should direct 
carrier flow normal to the sample thereby increasing the prob- 
ability of parallel (to electric field) spin-blocking (Figurel2tc)) 



as opposed to lateral (Figure |2|b)) spin-blocking (in plane of 
sample). Following previous researchers'^ '** ^^ we make the 
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FIG. 2. (Color online) (a) There are three vectors of importance: N, 
the vector normal to the sample surface, the applied field B, and dis- 
placement vector R. (b) R is lateral and (c) R is parallel to direction 
of current (N). In all, the two spins are denoted by their projections 
on the plane orthogonal to the applied field. 

assumption that the exchange and dipolar couplings are con- 
stants for each PP. A more general description of inhomoge- 
neous spin-spin couplings is an avenue for future work. 



III. PERCOLATION THEORY OF OMAR 

A key development in electronic transport theory in disor- 
dered media was the random resistor network model proposed 
by Miller and Abrahams over a half century ago.^^ In the sim- 
plest form where only spatial disorder is considered, the resis- 
tance between any two sites, / and j, is given by /?,■,■ = RQe^'''i''^ 
where r,y is their separation and I is the localization length of 
a carrier at a site. Miller and Abrahams were unsuccessful in 
determining the bulk resistance in their model; in fact it was 
not for over a decade that the resistance was calculated cor- 
rectly using a relatively new mathematical theory known as 
percolation theory! ^^"^^ ! The arguments of percolation theory 
dictate that bulk resistance is governed by a critical resistance 
(distance) Re (r^) which is the smallest resistance (separation) 
that permits an infinitely large chain of sites connected by re- 
sistors (separations) less than the critical one. The threshold 
length is constrained by the bonding criterion: 



Jo 



A%Nr^dr = Br 



(9) 



where A^ is the density of sites and Be tells how many bonds 
each site must connect to on average to be included in the 
percolating network; Be ~ 2.7 in three dimensions.^** Conduc- 
tion within r-percolation has been observed in organic semi- 
conductors in the regime of large inter-site separations and 
high temperatures where the influence of energy disorder is 
minimized. ^^ ^^ These conditions are assumed throughout this 
article. Despite the fact that transport in most organic semi- 
conductors may be determined by energy disorder alone,|^ 
our results elsewhere''' '^ and herein share many qualitative 
elements with OMAR experiments. This coiTespondence sug- 
gests that at least some aspects of OMAR may be independent 
of disorder type. 

For completeness, the spin-dependent percolation 



in many other OMAR theories, spin affects hopping transport 
through the Pauli exclusion principle; occupation of two po- 
larons at a single site (a bipolaron) is forbidden in the triplet 
state (T), but permitted in the singlet state (S) at an energy ex- 
pense U - the Coulomb interaction energy. For r-percolation 
alone, including t/ > provides a charge-blocking effect that 
diminishes spin-flip effects because bipolaron formation costs 
energy U regardless of the relative spin orientation. Since 
in reality energy disorder entails a spread in site energies 
(comparable to U) the charge-blocking effect is unphysical 
since there are frequently bipolaron energy states resonant 
with the polaron energy. This is shown schematically in 
Figure [3] So to avoid any unphysical charge-blocking effects, 
we make assumption that f/ = 0. A more sophisticated 
approach would examine the problem with energy disorder 
and/or variable range hopping^^, but is beyond the scope of 
this current work. 

A polaron with arbitrary spin is restricted from hopping to 
another polaron if the PP's spin state is T but may hop to an- 
other polaron (forming a bipolaron) if the result is an S state 
just as if to an unoccupied site. For a polaron looking to per- 
form a hop, the respective densities of these three types of sites 
are Nj, Ns, and A^o- Since the spins of forming PPs are ran- 
dom, spin statistics mandates that Ns — \Ni and Nt — |M , 
where A^i is the concentration of injected polarons. A^i is small 
compared to the total concentration of sites so that a polaron 
encountering a bipolaron is unlikely. 
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FIG. 3. Energy description of hopping between sites. Energy of a 
polaron/bipolaron is shown as a thick line. The dotted lines mark 
specific energy levels zero and —U. Solid circles represent actual 
carriers in space, (a) Polaron at L hops to unoccupied site R. Both 
sites are at equivalent energies, (b) Polaron at L hops to occupied site 
R which was at energy — (/. The Coulomb interaction energy of the 
bipolaron lifts the total energy back to zero. 



theory 
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is summarized here. In this theory, as also 



A PP that exists in the T state cannot form a bipolaron; 
this restriction effectively reduces the concentration of sites 
to N'eff = N — Nt since only these sites are accessible to a 
hopping polaron. When spin transitions are absent, the bond- 
ing criterion of Eq. d9 1 becomes A^ -^ ^'eff- ^^ne expects the 
resistance in the sample to increase since accessible sites are 
now typically further away from one another. Whether a PP 
combines into a bipolaron or disassociates as the polarons hop 
away from one another, the dynamics are integrally dependent 
on spin transitions since bipolaron formation is blocked for T 
states but unblocked for S states. The effects of spin transi- 



tions are now considered. Spin-blocked paths can be opened 
if the PP's spin function changes from T to S. A given polaron 
hopping to another site occupied by a polaron is successful 
1/4 of the time. The other 3/4 of the time the hopping attempt 
is rejected due to the PP being in the T state. At this point two 
things may occur: 

Bipolaron formation: The polaron waits until the next 
hop. If a spin transition (T-S) has occurred then the hop 
is successful. 

Disassociation: If the spin transition has not taken place 
and a more distant site is available, the polaron will hop 
to this further site. 

There is a trade-off between the two processes; if sites are 
well-packed or the hopping attempt frequency is very rapid, 
the polaron is likely to disassociate. If the opposite is true (di- 
lute and slow hopping frequen cy), th e polaron is prone to wait 
until the spin transition occurs .Elill Therefore spin transitions 
adjust the number of blocked sites from A^j- to (1 — pt^s)Nt 
where pr^s is the probability for the blockade to cease when 
the next hop is attempted. This probability may depend on 
time, spatial separation, applied field, and other system pa- 
rameters depending on what interaction initiates the spin tran- 
sitions. With the inclusion of spin transitions the bonding cri- 
terion becomes 
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4nNeffrdr = Be 



(10) 



with Neff =N — Nt+ pt^sNt- After defining the MR to be 
MR= {[Rc{B)-Rc{0)]/Rc{0)), we use Eq. ^ to write^ 



MR; 
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2^i y\ps{B)-ps{0))dy, 



(11) 



where we have made the substitution pr^s = i [1 ^ Ps] in 
terms of the easier to evaluate quantity ps^s = Ps and r| = 
Nt /N?^ This equation is valid only in the dilute limit (A^i <C 
A^) but this situation is believed to be the physically realized 
one. We have used a dimensionless length y — r/l. The 
renormalized critical length is jcj =ycg{l—NT/N)^^'^ where 
ycQ = (SBc/^Tif'Ny'^ is the critical length in the spinless 
problem, yc^ is nearly ycg because of the assumption of dilute 
polarons. In any case, we use the threshold as an independent 
variable and refer to it as simply yc hereon. Angular brack- 
ets denote an appropriate averaging which in our case is over 
hyperfine fields. 

The task of calculating the MR now reduces to finding the 
probabilities for singlets at time, t, that were initiated in the 
singlet state ps{B)- PPs do not experience a hop at exactly 
T/,; some hop earlier and some later and the probability that 
a hop has not occurred at time, t, is exp(—f/T/,).'2^ Therefore 
Ps{B) — ^ Jq ps{t)e^ ''dt where p is the two-spin density 
matrix. Spin-spin interactions affect the MR by altering the 
spin transition probabilities. 



We now consider the modifications required due to the new 
total Hamiltonian ,yf from Sec.|II]to obtain results for the MR. 
We must determine, from the total Hamiltonian, the singlet 
occupation probability at a time t given an initial singlet state 
(Ps). This was previously accomplished by using^^ 



Ps = |(5|expH^0|5)|' 



(12) 



Direct computation could follow after decomposing the expo- 
nentiated operator using the Euler formula. The key property 
of the Hamiltonian that allowed this process was that it was 
Zeeman-like. This is no longer the case with the spin-spin 
Hamiltonians now involving exchange interactions and dipo- 
lar terms. In fact no progress can be made using Ref. IT6l s 
technique since 



j{:Ho+:Hexc+J^d 



'dip)' a^ 



J^t i(:H,,,+^ji^)t 



(13) 



due to a lack of commutativity between Hamiltonians. If we 
wish to obtain analytic results, we are forced then to assume 
small spin-spin interactions and use perturbation theory. 



IV. EFFECT OF SPIN-SPIN INTERACTIONS ON MR IN 
THE FAST HOPPING LIMIT 

To ascertain analytic expressions in this section, we assume 
R\\N and that these vectors subtend an angle a from the ap- 
plied field which we take to be in the z direction. Then for 
small aeff, J, and D, but arbitrary cOq, the spin density ma- 
trix can be solved analytically using perturbation theory in the 
interaction representation. The non- intera cting case {J = 0, 
D — 0) has been examined in the pasF^'24J ^n^j -^q ^Iso solve 
it below as an instructional example. To the best of our knowl- 
edge, the interacting case has not been dealt with in this man- 
ner prior to the treatment here. 

In the interaction representation, the following operators 
are defined: 



M'{{t) 



oi^zt 



Me 



-i^zt 



(14) 



where Mi — Mhf + Mexc + Mdip is the perturbing Hamilto- 
nian. Similarly, the new spin density matrix is 



P*(0 



ol^Z' 



p(Oe 



-i.^ezt. 



(15) 



initially it can be shown that p*(0) = p(0). 

To simplify we consider the following: all terms odd in the 
hyperfine interaction vanish when the hyperfine interaction is 
assumed to be isotropic; additionally, it can be shown that 
the^rsf, second, and third order contributions from the spin- 
spin interactions do not affect the singlet occupation probabil- 
ity P5(r). The influence of exchange and dipolar interactions 
must be examined to fourth order perturbation theory. The hy- 
perfine interaction also contains a fourth order contribution. 

Several terms do not contribute and are neglected when we 
write out the perturbed density matrixl^^ 



p*{t)=p{Q)-fJ'^^ MO), j^i:f{t%jf:f{t")]dt"dt'+ 



+ :?4i^(r") + Hut'%.j^:f{t"') + :?4:,(r"') + Hun , j^^yin 



(16) 
dt""dt"'dt"dt'. 



A. Analytic solution in absence of exchange or dipolar 
interactions 

We begin by considering the solution in the absence of 
exchange or dipolar interactions. The hyperfine Hamilto- 
nian in the interaction representation is J^i*f{t) — U{J^tf^ + 



OqI 



2/5] sin 5M] 



5M 



Mf2)U^ with U = e'-^z' ^ [cos 

2!52sin^] since J^hfi commutes with •^hf^- We write 

p(0) = Ps where Ps is the singlet projection operator: 
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(17) 



in the 5, 7b, r+, T basis. 

The singlet part of the density matrix (5'|p*(f)|5') is ^s(t). 
Initializing in the singlet state requires that ps(0) = 1. The 
first order correction vanishes. After averaging over the hy- 
perfine fields for the ensemble of two carriers, including the 
second order term, Eq. ( [T6] ), yields 



(ps> = 1 - -Affi 



1+2 



sin^ coo? /2 

(coof/2)2 



(18) 



in agreement with the quantum mechanical calculation.E^ESl 

To find the MR, we substitute the singlet probability, Eq. 
( fTSJ l, into Eq. ([TTJ. After integrating over the exponential 
distribution of hopping times. 



MR: 



1 1 

3 y;:r^ 



COn 



CO, 



'0" 



1/t; 



y dy, 



(19) 



where r = {ihO-eff)^^ and zj^^ = voexp(— 2y). The integral 
in Eq. (fT9| can be computed when the hopping rate has 
an exponential dependence on the hopping distance T^y' = 
voexp(— 2^). The result is cumbersome but can be consid- 
erably simplified under the usual assumption that jc ^ 1 to 



MR: 



Yi'^h.^' 
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where r^ — (Xt-fle//) ' and x^ ' 
fields, MR,„,/ri = l/(12r2). 



voexp(— 23'c)- At saturating 



B. Analytic solution with exchange and dipolar interactions 

We now consider the addition of exchange and dipolar in- 
teractions. By using commutator algebra, the fourth order hy- 
perfine interaction can be separated from the spin-spin inter- 
actions which enter in quadratically with a quadratic hyperfine 



term. This implies that the spin-spin interactions would have 
no effect on the singlet probability if not for coexisting with 
the hyperfine interaction. This is understood by recognizing 
that the spin-spin Hamiltonians do not couple the singlet state 
with triplet states. 

The mathematics simplifies when no magnetic field is ap- 
plied. The second order correction can be obtained from Eq. 
( 18 1 in the appropriate limit. The time-averaged singlet prob- 



ability up to fourth order obtained from Eq. (16 1 is 



{ps{B^Q))^\ 



3 1 
4^ 



33 1 
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\> 
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(21) 



where the second term is the second order contribution and 
d — D/ueff, and j — J ja^ff. 

The effect of exchange and dipolar interactions is to lift 
the degeneracy of singlet and triplet states. In strongly ap- 
plied magnetic fields, the r+ and T states are far in energy 
from the S and 7b states due to the large Zeeman splitting. 
Consequently, only transitions between S and To states occur; 
the problem simplifies considerably since a smaller subspace 
of the original Hamiltonian warrants examination. The 2x2 
Hamiltonian matrix. 



-J (coi-co-2)/2 , 

(Wi-co|)/2 7 + D(l+3cos2a)/4y ' ^ ' 



is readily diagnolizable; the exchange splitting magnitude of 
S and 7b is CO^-t;, = aeff\d{l — 3cos^a) — 4y|/2. The time 
averaged singlet-to-singlet probability is 



{ps{B ^ -)) = 1 - 



11 3 1 1 1 «; 
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(23) 



^"ff 



To compare with Ref. l8 it is useful to define a quantity to 
measure the effects of angle variations: 



ln(l+cogT?)], (20) ®{B,a) 



{Ap{B,a)) _ {pT^s{B,a))-{pT^sm 
{ApiB,0)) {pT^siB,0)) - {pT^sm ' 



(24) 



The function is plotted for two values of the exchange con- 
stant in Figure |4] The number of maxima is contingent on the 
exchange coupling; the simulations of Wagemans et al. found 
that for J > — jD four maxima occur. '^ By examining our an- 
alytic expression for 0, the conditions of J that alter the num- 
ber of stationary points can be obtained from d&/da — 0. We 
find the analytic condition for the existence of the four max- 
ima. 



1 1 
-D<J < -D. 

2 - "4 



(25) 
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FIG. 4. (Color online) Angle dependence of singlet probability. Solid 
black: j = 0; dashed blue: j — —Id. For both, used d = 1 and r = 5. 
For ease of comparison with Ref . 1 8 t/, has been assumed constant. 



V. EFFECT OF SPIN-SPIN INTERACTIONS ON MR FOR 
ARBITRARY HOPPING RATES 

A. Calculation of the singlet probability 

As mentioned earlier, we are unable to use the analysis of 
Ref. [16] when spin-spin interactions are included. In the fast 
hopping we avoided the problem by using perturbation theory. 
However for arbitrary hopping times, we must use a different 
approach to finding the singlet probability. We begin again 



with Eq. ( 12 1; by inserting two complete sets of eigenstates of 
the full Hamiltonian the singlet portion of the density matrix 



Ps=Y. {S\&M~i'^t)\m){m\S){S\exp{Ly/^t)\m'){m'\S) 



The simulations of Ref. TW identified the existence of the 
lower bound, but here we report the full range where four 
maxima would be observed. The fourth order expression in a 
general field is available in closed form in the Appendix. We 
connect the occupation probabilities to the transport problem 
as described in Sectionjlll] For any magnetic field, the MR can 
be expressed in closed form, although not in a concise man- 
ner; it can be found in the Supplementary Information due to 
its length. 

At saturating fields. 
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(26) 

The full result is burdensome but without much error we can 
write down a simplified expression for the saturated MR: 



MR. 
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(27) 



where MRsaM = 0, ./ = 0)/r| = \/{\2rl). The angular depen- 
dence completely lies in the term Co|_j. causing the saturated 
MR to have different number of maxima depending on J and 
D just as with the case of 0. The angular dependence of the 
saturated MR can be examined analytically: 

MR,„,(90°) MR,„,(0°)_ 1 1 (coiji,(90°)-co2_ro(0°)) 



48 r4 
1 1 



^eff 



24,4^(4'^ + ^)' 



(28) 



where J can be seen to control whether the difference in re- 
sponses is positive or negative. Additionally, the angular re- 
sponse evaporates as hopping rate increases which leads us to 
the inference that spin-spin interactions and their concomitant 
angle-of-field dependence are unlikely to be observed experi- 
mentally in the fast hopping limit. 



£ e'(»„'-»m)' (5|,„) (,n\s) {S\m') {m'\S) 



= Y e'("'"'"""'^'|P?""'P 



where Ps = \S) {S\ is the S projection operator and its matrix 
representation is Eq. ( fTT] ). Since the desired quantity is ps — 
^ Io'Ps{f)^'dt, we first compute the time integral which 
can be readily accomplished: 



PS 



|pmm 



lAf 



(CO,, 



CO,, 



-l/v 



(29) 



This result is completely general for any Hamiltonian. This 
expression entails that only states with some singlet compo- 
nent contibute to ps- In this form, the spatial integral of Eq. 
( [TT] ) can be calculated yielding the result: 



ry'dy{psiB,y) = {'£\Pr'\ 
Jo \ .71,, 



yc Jcjj^Q , (C0„/-C0„)^ 

3 4 ^ v2exp(-4y,)^ 



Vc 



Li2(-..2 



1 . , («„ 

32^^^( 



-COm)^ 



1 

32 



(«„: 



(co,„/-co„,)2 
v2exp(~4y,.) 

2 



CO, 



Li3(— - 

vlexp{-4-yc)' 



(30) 



Li„(x) is the polylogarithm function of order n.l^The impor- 
tance of the polylogarithm terms dwindles as the hopping rate 
increases. In the case of slow hopping, such that vq -^ 0, only 
terms in Eq. ( 29 1 such that m' = m contribute to the MR which 

Dmm 1 2 \ 



(30ltO 'fiLml^S 



reduces the right-hand side of Eq. 

While the expression appears simpleTThe eigenvectors and 
eigenvalues of our Hamiltonian are complex enough that ana- 
lytic results are not possible. 

We first test this formalism against the results obtained in 
Ref. [16] The saturated MR, with spin-spin interactions absent, 
is plotted in Figure[5]and agrees with the prior formalism used 
in Ref. 16 

For the remainder of this article, the z direction is fixed to 
be parallel to N. We restrict R to lie parallel to N except where 
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FIG. 5. (Color online) MRjo/ as a function of hopping rate with no 
spin-spin interactions. Black solid circles: y^ = 5; red triangles: Vc = 
7. Solid lines are analytic result, Eq. [27] from 4th order perturbation 
theory. Inset is same as main except on log-linear axes. Averaging 
required 10^ configurations of the hyperfine fields. Results agree 
with alternate method of Ref . J^ 
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FIG. 6. (Color online) Saturated MR versus exchange constant for 
different values of the hopping rate; solid black: vo/a^ff = 0, dotted 
blue: vo/agff = 2 x 10 , dashed orange: vo/a^ff = 10 . For all, 
yc = 5. Inset: MR in the slow hopping limit for several different 
values of the exchange constant: solid red: J = 0.4, dotted green: 
y = 1, dashed purple: J = 3. For all, yc = 5. 



we allow R to vary over a range of angles. We choose to iso- 
late the effects of the exchange and dipolar interactions by ex- 
amining them independently in this and the next subsection. 



In the final subsection V D both interactions are analyzed in 
tandem. 



B. Exchange interaction 

Figure l6] plots the saturated MR versus exchange constant 
for several different hopping rates. As expected from the non- 
interacting picture, the overall MR decreases with hopping 
rate. When we concentrate on the zero hopping limit (solid 
black line in Figure l6]l, it is noteworthy that at small and large 
exchange couplings the results agree well with the Figure [5] 
for small and large hopping rates. 

We interpret this as the exchange interaction acting much 
like hopping; instead of the carrier hopping site to site, the 
spin of the carriers swaps between occupied sites. This heuris- 
tic understanding is further supported by the calculation of 
{ps{B — > oo)), wherein we find exchange can be accounted 
for exactly by the substitution 1/x/, — )• l/x/, + 27. If J/a^ff is 
small, the hyperfine fields fully mix spin before the spins swap 
and the case is no different than for slow hopping and no ex- 
change. If J/oeff is large, spins swap back and forth faster 
than hyperfine fields can mix which leads to a motional nar- 
rowing effect much like what occuiTed when the hopping was 
very rapid (but hopping is effectively zero here). This effect 
has been used in doped inorganic semiconductors to explain 
certain long spin lifetimes due to motional narrowing by ex- 
change interactions.^*^ ^^ A non-zero hopping rate suppresses 
the overall saturated MR as discussed in Ref. 16 

The large MR enhancement prior to J/aeff = 1 in Figure 
[6] is explainable in the following way: as J increases from 
zero, both ps{B — )■ o°) and ps{B = 0) increase sharply when 
/ ^ Ueff since the amount of spin mixing will be significantly 



reduced when the exchange interaction surpasses the hyper- 
fine interaction. The enhancement is due to the fact that both 
functions do not increase in the same way; ps{B = 0) in- 
creases more slowly with J than ps{B — > oo). As 7 increases 
far past flg//, the zero and infinite field probabilities become 
nearer in value (unity) since, as in the case of faster hopping, 
the hyperfine fields produce less spin mixing in both cases. 
The reason for the lag in ps{B = 0) is readily apparent when 



considering the pertinent energy levels and Eq. (29 1. At zero 



field, all four electronic energy levels are roughly degenerate 
when J ^ aeff. In a large field, two energy levels are already 
well separated and lead to negligible spin mixing. Therefore 
it takes a larger J to suppress spin mixing in zero field than it 
does in large field. 



Eq. 29 shows how energy level placement affects spin mix- 
ing; energy levels well separated from the S level do not factor 
into spin mixing (and pg increases toward unity). In the ab- 
sence of exchange, increasing the applied field separated out 
the r+ and T states from the Tq and S states which is why 
MR>0 (since ps{B) > psi^))- Inclusion of exchange allows 
for an USEE: when coq ~ 2\J\, the T+ or T levels are near 
the S or 7b levels, which increases the amount of spin mix- 
ing such that ps{B) < psiO) (which leads to MR<0). MR 
is shown for three different exchange constants in Figure l6|s 
inset. For finite yet small J (black solid line), the USEE de- 
velops and becomes more pronounced as J further increases 
(dotted green and dashed purple lines). The exchange inter- 
action may serve as the explanation of the USEE measured 
in unipolar organic devices.''* Increasing the hopping rate de- 
grades the large negative MR produced by the large exchange 
interaction (not shown). 



C. Dipolar interaction 

In this section J is set to zero and the dipolar interaction 
constant D is turned on. We focus here solely on the slow 
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hopping limit. Figure IT] depicts the dependence of saturated 
MR versus D/ueff on the angle-of-field, a. The difference is 
dramatic for D/ueff ^ 1. At angles ai^ound jCOS"^(— j) w 
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FIG. 7. (Color online) Saturated MR plotted at a = 0° (solid black), 



45° (dotted blue), a = ^cos '(- 



55° (dash-dotted red). 



and a = 90° (dashed orange). Hopping rate and exchange are zero. 
Non-zero hopping rate reduces overall MR^a/. Inset: same as main 
except polaron pair displacement vector R is not fixed but is averaged 
over an angular range within 10° of the z- 

55° the saturated MR goes strongly negative. At the angle 
5 cos^ ^ (~ 3 ) '^he dipolar splitting vanishes in the high field re- 
gion; a unique situation is created in which the zero-field split- 
ting of the energy levels is larger than the high field splitting 
between the To and S levels. This ensures more hyperfine mix- 
ing and less resistance at high fields than at zero field which 
lead to negative MR. Such large angular dependence has not 
been observed. This may be explainable by the inset of Fig. IT] 
where it can be seen that the special angle loses its uniqueness 
when the dipolar coupling is allowed to vary within 10° of 
N. Therefore it is unlikely that the special angle could be ob- 
served in experiments since the displacement vectors between 
PPs are not fixed but vary over a range of angles. 

The angular dependence is plotted as dashed red curves in 
Figure [8] from a = 0° to 180°. The concavity of the blue 
(solid) curves is set by the sign of the exchange interaction; 
hence, if the angle dependence is due to spin-spin interactions, 
angle measurements of OMAR could indicate the sign of the 
exchange interaction if the direction of the dipolar anisotropy 
is established. In the two insets, the strong angular depen- 
dence due to the special angle, jCos^'(— j), is seen to di- 
minish as the PP orientations are made to be less anisotropic. 
In the limit that the PPs orientations are completely isotropic 
(not shown), the angular dependence vanishes as expected. 



D. Exchange and dipolar Interactions 



In the regime opposite to the one considered in Section IV 
(i.e. slow hopping/large spin-spin coupling), we numerically 
determine that the same conditions on J and D found in Sec- 



maxima in the angle (0° to 360°) versus MR.sa, graph as shown 
in the main portion of Figure [8] These maxima continue to 
occur since the high field dipolar splitting varies with angle 
which in turn varies the amount of S — Tq mixing; however the 
modulation now depends on / and D: a = j cos^^ (— 5 — |g) 
and is not pinned near 55° as was the case without exchange. 



Pi 




tion IV Eq. ( 25 1, are still valid and produce either two or four 



a 



FIG. 8. (Color online) Main: sohd blue: rf = 0.1, j = 0.2; dashed 
red: d = 0.1; 7 = 0. Inset left: solid blue: d = 0.2, j = 0.2; dashed 
red: d = 2, j = 0. Inset right: same as left inset except polaron 
pair displacement has an angular spread around z of 35°. 4 x 10 
hyperfine configurations were used in the averaging procedure. 

An important conclusion to draw is that if the angle de- 
pendence of MR is due to dipolar coupling, the exchange in- 
teraction must be present and non-negligible since zero ex- 
change yields the wrong number of maxima. Neither can the 
exchange term be much greater than the dipolar term since 
it will suppress MR in that case (not shown). So the asser- 
tion that the exchange and dipolar couplings affecting MR are 
comparable seems necess ary d espite the drastic difference in 
their spatial dependencies.! P" ^ As already stated, a dipolar in- 
teraction initiated angle dependence is contingent on there be- 
ing a PP orientation anisotropy. 



VI. CONCLUSIONS 

Recent experiments have measured an angle-of-field depen- 
dent magnetoresistance. One possibility for this effect is an 
anisotropic dipolar coupling between hopping polarons. To 
evaluate the role of dipolar coupling in these systems, we 
have calculated magnetoresistance in positionally disordered 
unipolar organic systems. In accord with other authors,'** the 
magnetoresistance response to the dipolar interaction must 
be examined in conjunction with exchange coupling to agree 
with experiments. Sorting out the cause of the magnetore- 
sistance dependence on angle-of-field requires experimenta- 
tion in more organic systems that may allow for variations in 
the spin-spin couplings. Observation of four maxima would 
be unequivocal evidence for the dipolar interaction being the 
source of magnetoresistance variation due to angle. Another 
possibility for the angle dependence is the existence of an 
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anisotropic hyperfine interaction. The theory described here 
can be modified to account for anisotropic hyperfine distribu- 
tions and is the subject of future work. 

The theory outHned above is for unipolar organic devices. 
Many experiments are believed to be in a bipolar current 
regime. We previously argued^^ that our theory is still qual- 
itatively useful in describing bipolar MR with the relation 

MR„„,po/„,- ^ —A- MRjip,jiar where A, is some proportionality 
factor. This is identical to saying MRj„„p„/„r ^ ^ ^Cdipoiar- 

In summary, we have calculated the effects of spin-spin 
interactions on magnetoresistance in positionally disordered 
organic semiconductors. The magnetoresistance shows an 
angle-of-field dependence which agrees with experiments for 
certain ratios of the exchange and dipolar parameters. Future 
experimentation is needed to decide whether the angle depen- 
dence is due to the dipolar coupling or some other anisotropic 
mechanism. In addition to the utility of these calculations to 
the fledgling field of organic spintronics, we remark that these 
calculations of singlet and triplet probabilities including both 
the exchange and dipolar interactions may find application 



in the study of magnetic field effects on chemical reactions 
and more sp ecific ally their manifestation in avian navigation 
mechanisms .EIH2| 
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Appendix A: 4"* Order Singlet Probability 

The fourth order singlet probability is composed of a hy- 

(4) (4) 

perfine only term (p\ ) and a mixed term (p2 ) of hyperfine 
and spin-spin couplings: p^ = Pi - 



P2 where 



and 



(4)\ 



24/j^ + 90/zV2 + 145/!^/ + 136/iV*^ + 33r** 



(Al) 



ipi') - 



:2\ A 



h^ {nd^ + 32dj + 128/) + 3h'^ (1 lc/2 + 32dj + 128/) r^ +A%h^ {d^+4-dj - 8/) r 
128r^ (/z2 + r2) L 

48 {d^ + 8f) r^ + ndh^ {h'^{d + 8j) + 3h^{d + 8j)r^ +48;>'*) cos(2a) +9<i2/i4 [h^ + 3r^) cos(4a) 



(A2) 



with h = (Oo/aeff. MR can be calculated from these equations 
but due to their length, they are located in the Supplementary 
Information. 
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